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Abstract 

Several galaxies, such as dwarfs and Irregulars as well 
as outer galactic clouds have low metallicity. At low metal- 
licities a reduction in the amount of dust and heavy el- 
ements plays a significant role on the chemistry as well 
as the heating and cooling of the gas in the molecular re- 
gions, called as Photon Dominated Regions (PDRs). We 
present here the effects of reduced metallicity in PDRs 
and study the important PDR cooling lines ([CII], [CI] 
and CO). Moreover many observational evidences suggest 
that molecular clouds are clumpy. We model the molec- 
ular emission from galaxies incorporating a mass spec- 
trum of clumps. We also compare our results with the 
semi-analytical results obtained by Bolatto et al. (1999). 
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1. Introduction 

Photon Dominated Regions are predominantly molecular 
and atomic regions where the physical and chemical pro- 
cesses are dominated by UV radiation (cf. HoUenbach & 
Tielens, 1997, 1999). The molecular clouds in the vicinity 
of the newly formed stars is heated by the FUV photons 
in the energy range from about 6 to 13.6eV. These clouds 
cool through the atomic and molecular spectral lines, 
such as [Cll]l58nm, [OI]63,146Aim, [Cl]609,370nm and 
the milli-metric and sub-mm CO rotational lines. Plane 
parallel and spherical models of PDRs have been con- 
structed to understand the physical and chemical char- 
acteristics of these regions (e.g., Kaufman et al., 1999; 
Koster et al., 1994; Le Bourlot et al., 1993; Sternberg & 
Dalgarno 1995; Storzer et al., 1996; Tielens & HoUenbach 
1985). However there are other important factors which 
affect the UV absorption and scattering as well as the 
basic heating and cooling processes in PDRs. 

Several galaxies such as Dwarf galaxies. Irregular galax- 
ies, the Large and Small Magellenic Clouds have low 



metallicity (cf. Wilson 1995). A radial gradient in metal- 
licity of molecular clouds is found in the Milky Way and 
several other nearby galaxies (e.g., Arimoto et al., 1996). 
These low metallicity systems have much higher [CII] /CO 
and [CI] /CO line ratios as compared to the galactic value 
(e.g. Bolatto et al., 2000; Madden 2000). This suggests 
that the effects of metallicity should be considered while 
interpreting the molecular and atomic spectral line ob- 
servations of these sources. Since the important surface 
coolant of the PDRs, the [CII]158/Ltm emission, is strongly 
affected by the metallicity factor, we study the effects of 
metallicity in PDRs. In addition, low metallicity PDR 
models would also help us to understand the star forming 
regions in Dwarf galaxies which resemble the primordial 
galaxies. 

Additionally, observations of edge-on PDR have sug- 
gested that the molecular clouds are clumpy, and the UV 
radiation can penetrate deep inside the clouds (cf. Stutzki 
et al., 1988; Boisse 1990). This suggests that a simple 
single component model may not be sufficient to explain 
the observed features. We also include a mass spectrum 
of clumps to understand the cooling lines of PDRs from 
low metallicity galaxies. 

2. Metallicity effects 

Low metallicity systems have a lower, dust to gas ra- 
tio and heavy elemental abundances as compared to the 
galactic ISM. This reduction in the amount of dust grains 
affects, the absorption of UV radiation, heating of gas by 
photo electric emission from dust, formation of H2 on the 
dust grains and the cooling of the gas through atomic 
and molecular lines. In addition, the chemistry is also af- 
fected by the reduction in dust and heavy elements (van 
Dishoeck & Black 1988; Lequeux et al., 1994). We use 
the self-consistent spherical PDR model of Storzer et al. 
(1996) to study the metalhcity effects. We scale the dust 
dependent parameters and the abundance of heavy ele- 
ments with the metallicity factor, Z, in our PDR calcula- 
tions. 
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The variation of the temperature at the surface of 
the PDR clumps with the incident UV field is plotted 
in figure n for a clump of mass, M=1M0 and density, 
n=10'^cm~'^. The UV field, x, is expressed in units of 
mean UV field of Draine (1978). It is seen from the figure 
that at high UV fields the temperature is proportional to 
metallicity, whereas at low UV fields there is no significant 
change in the surface temperature. 
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Figure 1. The temperature at the surface of the PDR 
clump of mass M — 1 Mq and density n = l{ficm~^ 
is plotted against the UV field for different metallicities. 



This correlation can be understood analytically by 
balancing the dominant heating and cooling processes. In 
the case of PDRs exposed to high UV fields grain photo- 
electric emission (PE) dominates the heating. The photo 
electric heating rate given by Bakes & Tielens (1994) is 



TpE = lQ'^\nZ 



3 X 



1 + 2 X 10 



(1) 



where Z is the metallicity factor. 

The cooling is dominated by fine structure [01] emis- 
sion, gas-grain collisions and fine-structure [CII] emission. 
The cooling rate can be written as Lij — niAijhvijl3{Tij) 
where (3 is the escape probability, Aij is the transition 
probability and number of atoms at level i and huij is 
the corresponding frequency of radiation. 

Following HoUenbach & McKee (1979), for the popu- 
lation of the first level of the 01 atom, the cooling rate of 
[OI]63/xm can be written as. 
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Figure 2. The equilibrium temperature for the balance be- 
tween grain photo-electric heating rate and [01] cool- 
ing rate is plotted against the metallicity, for a density, 
and UV field, x = 10^. 
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8.45 X lQ-^'^l3{T)nHZ 
l + 5e228/T + i/3e-98/T- ^ ) 

Although eqns. Q and El show a similar dependence 
with metallicity, the grain heating rate has an additional 
dependence on metallicity through the charge state of the 

grains. The charge state of the grains is expressed as — 
which depends on the availability of electrons for recom- 
bination. The main source of electrons at high UV fields is 
the ionisation of CI. Since at the surface almost all of the 
CI is ionised, the electron density rig ^ nc oc nc^z=iZ. 
With this assumption the heating rate, Vpe/Z decreases 
with decreasing Z whereas the cooling rate, LqI/Z re- 
mains constant. The equilibrium temperature obtained 
by solving eqns. and El decreases as the metallicity de- 
creases (cf. figure \^ as seen in our PDR calculations. 

3. Clumpy PDR model 

In our clumpy model, the molecular cloud is modelled as 
being composed of many spherical clumps of mass spec- 
trum of the form, = AM'". We use a = 1.8 (e.g. 
Kramer et al., 1998). We assume that the turbulent ve- 
locity dispersion of the cloud is larger than the intrinsic 
line width of each clump. Thus the clumps do not interact 
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radiatively and the total intensity of a spectral line from 
the cloud can be written as, 

It - / I{M) — ft,c fM,bdM 

JM^i,, cloud 

The beam filling factor of each clump is fM,b{— ^Ai/^b) 
where Hm is the solid angle of the clump of mass M and 
fif, is the beam solid angle. The fraction of the clumps 
within the beam is given by /f,,c — ^t/^c- 

By scaling the mass with Mmax {x — M / M„iax) 



A « M,. 



a)M, 



Mrr 



(2-a) 

max 



/M„ 



IxX-^f. 



M,Mrr 



dx 



where is the beam filling factor, xq — x{Mjnin) and 
fx = R^{x)/R^ig^j^- Although the total intensity depends 
on the scaling constant and the beam filling factors, the 
line ratio between any two spectral lines depends only on 
Mmtn/Mmax and a. 

4. Clumpy clouds at low metallicity 

The observations of [CII] and [CI] emission from low metal 
galaxies have been modelled by Bolatto et al. (1999), as- 
suming that the size of the CII region scales inversely 
with metallicity. It is also assumed that the size of the CI 
region remains between two limiting scenarios of, an in- 
verse dependence and no variation with metallicity. Our 
spherical PDR model calculations based on the model of 
Storzer et al. (1996) show that the size of the CII layer 
is indeed inversely proportional to the metallicity factor 
Z. However the size of the CI region shows very weak de- 
pendence at low Z and a roughly inverse dependence at 
high Z for a typical spherical clump of density, lO'^cm"'^. 

The intensity ratio I[cii]i58tJ.m/ Ico{i-o), observed in 
many nearby low metal galaxies, show a power law de- 
pendence with metallicity. This dependence has been pre- 
dicted using a semi-analytical clumpy model by Bolatto 
et al. (1999). We use our clumpy model explained in sec- 
tion O to study the metallicity dependence of this line 
ratio. The observed ratios of nearby galaxies can be well 
represented by a clumpy model of density 10^ cm~^, ex- 
posed to a UV field of x = 10^ as shown in figure 01 
The higher observed ratio for the 30 Doradus region can 
be explained by a similar clumpy model, but exposed to 
an UV field of x = 10^. These results compare well with 
the results shown by the semi-analytical model of Bolatto 
et al. (1999). This trend suggests that at low metallici- 
ties CII is a tracer of molecular hydrogen rather than CO. 
However the observed variation of [CI] /CO line ratio with 
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Figures. The intensity ratio /[c//]i58MmAco(i-o) 
plotted against the metallicity. The maximum and mini- 
mum masses of the clump spectrum are l(fi and 10^^ Mq. 
The density at the surface of the clumps are l(fcm^^. 
The observed ratios of nearby galaxies are plotted as filled 
squares. The model result plotted as solid line compares 
well with that of the trend shown in figure 5 of Bolatto et 
al. (1999). 

metallicity is steeper than the model prediction. This is 
most likely due to large [CI] line intensities predicted by 
the PDR models and requires further investigation. 
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